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Abstract 

The  reaction  probability  per  collision,  Pc  ,  for  1 -methyl cyclobutene 
activated  at  a  hot  seasoned  fused  silica  surface  has  been  measured  under 
single  collision  conditions  over  the  reactor  temperature  range  Tr  = 

480  K  -  800  K.  The  initial  thermal  vibrational  energy  population  distri¬ 
bution  of  the  cyclobutene  molecules  was  varied  from  Tc  =  273  K  -  570  K. 

The  reaction  is  isomerization  to  isoprene.  Mixed  samples  of  1-methyl- 
cyclobutene  and  cyclobutene  (which  had  been  studied  earlier)  were  used 
for  internal  comparison  under  the  same  conditions.  These  experiments 
provide  a  test  of  the  relative  suitability  of  various  analytical  forms 
for  the  collisional  transition  probability  matrix  P  .  Stochastic  calcu¬ 
lations  with  a  Gaussian  form  provide  the  best  overall  fit  to  the  data. 

The  calculated  average  amount  of  energy  (<AE'>r  )  transferred  from  the 

fc0 

hot  molecules  in  a  vibrational  down  transition  from  the  reaction  threshold 
energy  level.  Eg  ,  declined  from  7220  cm’^  to  3890  cm"^  with  increase  in 
surface  temperature  from  600  K  to  800  K.  The  experimental  collisional 
efficiency,  8^  ,  declined  from  0.39  to  0.035  over  the  combination  temper¬ 
ature  range  Tr,Tc  =  600,500  to  800,293.  Strong  collider  behavior  was 
observed  with  both  1 -methyl  cyclobutene  and  cyclobutene  for  T  less  than 
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Introduction 

Gas-surface  vibrational  energy  accommodation  is  being  intensively  studied, 

1-5 

both  experimentally  and  theoretically.  The  variable  encounter  method  (VEM) 
provides  a  simple  technique  for  the  study  of  colli  si onal  transfer  of  vibra¬ 
tional  energy  between  gas  molecules  and  a  surface  at  levels  of  excitation 

corresponding  to  homogeneous  unimolecular  reaction.  By  this  technique,  cold 
substrate  molecules  that  are  initially  equilibrated  at  low  temperature  Tc  experi¬ 
ence  a  known,  and  experimentally  variable  number  of  sequential  collisions,  m  , 
with  a  hot  reactor  surface  at  temperature  before  leaving  the  reactor  and 
re-equilibrating  to  their  initial  low  temperature.  Values  of  m  between  2  and  30 
have  been  conventionally  used.  Gas  molecules  describe  a  random  walk  along  an 
energy  coordinate  until  they  reach  an  absorbing  level,  the  critical  reaction 
threshold  EQ,  characteristic  of  the  homogeneous  gas  reaction.  Such  relaxation 
of  non-equilibrium  vibrational  energy  distribution  was  described  theoretically 
a  number  of  years  ago  by  Rubin  and  Shuler,10  and  by  Kim11  and  Widom.12 

Recently,  single  collision  measurements  (m =  1 )  by  the  VEM  technique  were 

described  for  the  cyclobutene  (CB)  system.12’1^  Initial  vibrational  energy  of 
14 

CB  was  varied.  In  the  present  communication,  we  have  applied  the  single  colli¬ 
sion  condition  to  the  study  of  1 -methylcyclobutene  (1-MCB)  system  for  comparison 
with  that  of  CB.  The  surface  is  a  "seasoned"  fused  quartz  finger.  Although  such 
a  surface  is  not  well  defined,  it  is  the  conventional  experimental  surface  of 
thermal  kinetics. 

The  vibrational  energy  population  vector  of  molecules  after  one  collision, 

N1  ,  is  given  by  N1  =  PNc  eq,  where  P  is  a  collisional  transition  probability 
matrix,  and  Nceq  is  the  initial  vibrational  energy  population  vector  that  corres¬ 
ponds  to  the  thermal  Boltzmann  distribution  at  the  low  temperature  of  the  gas 
reservoir  wall.  Since  no  comprehensive  theory  of  gas-surface  collision  inter¬ 
action  exists,  at  least  for  complex  molecules  on  these  surfaces,  down-jump 


transition  probability  elements  of  the  P  matrix  were  constructed  according 
to  various  assumed  mathematical  models;  these  have  some  plausible  connections 
with  physical  reality.  Up-jump  transition  probabilities  were 
constructed  from  the  down-jump  transitions  with  use  of  the  conditions  of 
detailed  balance  and  completeness. 

The  ring-opening  unimolecular  isomerization  of  1-MCB  to  isoprene  has  a 
low  Eq  { 34.2  kcal  mol"^);  that  for  CB  is^5  32.4  kcal  mole”^.  The  effect  of 
the  methyl  substituent  on  the  collisions!  reaction  probability  provides  an 
interesting  comparison.  Previously,  a  similar  comparison  was  made  for  the 
collisional  relaxation  of  vibrational  energy  transients  in  the  methyl  cyclopropane 

g 

and  cyclopropane  systems  using  a  VEM(m  >  1)  technique. 
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Experimental 

The  1-MCB  substrate  sample  was  a  mixture  containing  14%  of  3-MCB  as 
impurity;  it  was  not  removed  since  it  caused  no  complication  and  allowed  a 
rough  concurrent  measurement  for  its  own  isomerization.  A  similar  amount  of 
cyclobutene  was  added  as  an  internal  standard.  The  reaction  system  apparatus 
used  here  was  the  same  as  that  for  previous  single  collision  experiments.13,14 
The  reaction  vessel  consisted  of  a  3-5.  spherical  pyrex  reservoir  flask  that 
was  provided  with  an  internally  heated,  central  fused-quartz  finger.  The  flask 
and  the  finger  were  heated  independently.  Two  experimental  series  were 
used  to  measure  the  isomerization  rate  constants:  a)  the  surface  temperature 
of  the  reactor  finger  was  varied  from  600  K  to  800  K,  while  the  wall  temperature 
of  the  flask  reservoir  was  varied  from  273  K  to  550  K.  b)  the  reactor  tempera¬ 
ture  was  varied  from  480  K  to  800  K  at  constant  wall  temperature  of  273  K. 
Temperature  deviation  of  the  finger  surface  was  a  maximum  of  t  5°K  at  800  K, 
and  that  of  the  flask  wall  was  t  10°K  at  550  K.  Before  kinetic  measurements 
were  made,  the  reactive  surface  of  fused  quartz  was  "seasoned"  at  the  highest 
temperature  and  the  seasoning  was  maintained  by  exposure  for  a  few  hours  to 

_3 

the  mixture  gas  at  a  pressure  of  'v 5x10  torr  prior  to  each  run. 

The  system  was  run  in  both  static  and  flow  modes  in  the  pressure  region 
-4  -3  -4 

between  10  and  10  torr,  usually  'v  2x10  torr.  For  the  flow  mode,  typical 

residence  times  in  the  reactor  were  10  to  30  sec.  Duration  of  a  flow  run  was 

13 

several  minutes.  Detailed  description  of  both  modes  was  given  earlier. 

Pressure  measurements  were  made  with  an  MKS  146H  capacitance  manometer. 

Product  analysis  was  performed  by  gas  liquid  phase  chromatography  on  a 
5  ft  x  3/16  in.  squalane  column  on  Chromsorb  P  at  0°C  with  FID  detection. 
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Results 


Observed  rate  constants  for  1-MCB  isomerization  to  isoprene  were  calcu¬ 
lated  from  the  product  yield.  The  total  observed  rate  constant  was  a 
sum  of  reaction  due  to  heating  at  the  reaction  finger,  k^  ,  and  at  the 
reservoir  wall,  kc  ,  given  by 

“t'VV  ■  kr(VTc>  *  kc(Tc> 

Independent  measurements  of  kc(T  )  were  made  under  the 
experimental  condition,  Tr  =  Tc  »  with  use  of  a  minor  area  correction  for  the 
relative  surface  area  of  the  reactor  finger  and  the  reservoir  wall  (1:13). 

Both  rate  constants  kt  and  k£  are  shown  in  Fig.  1.  Experimental  values  having 
a  cross  mark  in  Fig.  1  were  rejected,  since  correction  for  k  amounted  to  =  50% 
of  the  total  rate  k..  The  reaction  probability  per  collision,  P  (T  ,T  ),  was 

t  L  I  L 

calculated  from  the  finger  surface  reaction  constant  kr(Tr,Tc) , 

Pc(Tr*Tc)  =  kr(Tr,Tc)/(A(8kTc/nm)i/4V)  , 

where  A  is  the  reactor  finger  surface  area,  k  is  the  Boltzmann  constant, 
m  is  the  molecular  weight,  and  V  is  the  volume  of  the  reservoir  flask.  The 
values  of  the  reaction  probability  are  plotted  as  a  function  of  Tc  in  Fig.  2 
and  are  listed  in  Table  1.  The  experimental  uncertainty  in  P£  is  estimated 
to  be  20%.  Also  listed  is  the  collisional  efficiency,  B-j  ,  defined  as 
$1  =  pc{Tr»Tc)/pc  (Tr);  these  quantities  are  the  analog  of  the  conventional 

C  r 

homogeneous  thermal  collisional  efficiency  factor,  B  .  The  values  of  Pc  (Tr) 
are  calculated  from  the  Boltzmann  distribution  population  vector  characteristic 


Jfl 


The  mathematical  forms  used  here  for  the  transition  probability  distri¬ 
bution  of  down-transition  energy  jumps  AE  are  exponential  (E),  Gaussian  (G), 
a  Boltzmann  weighted  exponential  (BE),  and  Gaussian  (BG)  functions,  given  by 
Eqs.  ( 1 ) - ( 4 ) : 

pE(AE)  =  C1  exp  (-AE/<AE>)  (1) 

pG(AE)  =  C2  exp  (-(AE-AEmp)2/2o2)  (2) 

P^j(AE)  =  C3BiPE(AE)  (3) 

pSj(AE)  =  C4BiP6(AE);  Bi  =  g  .exp(-E./RTr)  (4) 

Here,  <AE>,  AEmp,  and  o  are  parameters  of  the  models,  the  C^s  are 

normalization  constants;  and  B.  is  a  normalized  Boltzmann  distribution 

characteristic  of  the  finger  surface  temperature,  ,  where  gi  is  the 

density  of  internal  states  at  energy  level  E.;  o  is  set  equal  to  0.7  AE 

l  mp 

for  the  G  and  BG  models.  Transition  probability  values  below  E  =  0 

in  a  down  transition  from  the  energy  level  i  were  added  to  the  element  for 

elastic  collisional  transition  probability,  p^  .  The  E  and  G  models  have 

been  termed  "flat"G’^  since  the  distributions  given  by  them  are  independent 

of  the  initial  energy  level.  The  microscopic  rate  constants  at  each  reactive 

energy  level  were  calculated  by  RRKM  theory  with  use  of  the  vibrational  fre- 
.  15 

quency  assignment  of  Elliott  and  Frey;  these  calculations  were  made  with 
an  energy  grain  size  of  100  cm’1  that  was  also  used  to  specify  the  transition 
probability  matrix.  Isomerization  of  excited  molecules  occurs  during  the 
average  (collisionless)  flight  time  between  the  reactor  surface  and  the  reser¬ 
voir  wall  (9  cm).  An  average  fraction,  f^  ,  of  molecules  excited  above  EQ 
that  decompose  during  the  flight  time  was  calculated  in  order  to  estimate 
the  dependence  of  the  observed  rate  constant  on  the  reaction  probability. 

The  most  favorable  experimental 
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situation  is  obtained  for  =  1,  corresponding  to  k/k^  -*■  0  in  a  homogeneous 
system;  in  this  case,  which  holds  for  CB,  reaction  is  governed  only  by  the  colli- 
sionally  activated  population  distribution,  and  the  need  for  RRKM  calculations  of 
k(E),  the  specific  reaction  probability  at  internal  energy  E,  with  a  postulated 
activated  complex  structure,  is  removed.  In  the  present  case,  f^  ^  0.6-0.65,  so 
that  the  accuracy  of  the  measurements  was  only  moderately  dependent  on  the  accu¬ 
racy  of  the  RRKM  calculational  details.  Clearly,  as  f^  0,  corresponding  to 
k/k^  1  in  the  p  -*■  <*>  homogeneous  case,  no  reliable  information  about  energy 
transfer  can  be  gained. 

Results  of  the  least  squares  fitting  to  the  experimental  P£  curve  for  each 

reaction  surface  temperature  (with  use  of  <AE>,  or  AEmp,  as  the  parameter  of  fit) 

are  shown  in  Fig.  3  with  use  of  four  models.  The  corresponding  parameters,  the 

average  down-transition  energy  <AE>  for  E  and  BE,  and  the  most  probable  energy 

AEmp  for  G  and  BG,  are  listed  in  Table  2.  We  have  chosen  to  exhibit,  and  use, 

these  average  values  rather  than  to  enumerate  a  plethora  of  best-fit  values  for 

each  value  of  T£  for  a  given  Tp.  Obviously,  use  of  an  average  constant  parameter 

for  the  whole  Tc  range  cannot  give  as  good  fit  at  each  Tc  as  would  an  optimized 

value.  The  calculated  value  of  <AE'>F  is  the  average  amount  of  energy  trans- 

t0 

ferred  in  a  down-transition  from  the  threshold  energy  level  Eg;  <AE‘>  is  independent 
of  Eg  for  E  and  G,  but  varies  with  initial  level  for  BE  and  BG.  Also  listed  is 
<AE+>av  =  >Eq'  Eo]’  the  avera9e  amo'jnt  of  energy  transferred  in  an  up 

transition  to  the  levels  above  Eg.  The  average  vibrational  energy  of  1-MCB 
molecules,  E  ,  for  the  thermal  Boltzmann  distribution,  is  given  for  the  temper¬ 
atures  of  interest  in  Table  3.  Another  important  quantity  is  the  average 
energy  transfer,  AEflv,  defined  as  AEflv  =  Ef  -  Ec,  where  Ec  and  Ef  are  the  aver¬ 
age  energies  of  molecules  before  and  after  single  collision;  AEav  is  listed 
from  results  of  E,  G,  and  BE  calculation  in  Table  4.  The  vibrational  energy 

accommodation  coefficient,  a  ,  is  defined  as  a  =  AE  /(E  -E  ),  where  E  is  the 

flv  r  c  r 

average  energy  E  corresponding  to  T  =  Tr  and  is  given  in  Table  4. 


Single  collision  experiments  permit  a  more  refined  test  of  the  correct 

form  of  P  than  do  m  >  1  experiments.  In  particular,  an  even  more  stringent 

criterion  is  realized  here  by  changing  the  initial  population  vector  N  .  The 

features  of  the  experimental  Pc  curves  for  all  Tr  in  Fig.  2  are  similar  to 

those  for  CB  and  also  exhibit  a  foot  for  Tc  <  400  K.  The  best  fits  at  higher 

temperatures  were  obtained  with  the  E  and  G  models.  The  BE  and  BG  models 

fail  to  agree  with  the  data  curve.  At  lower  temperature  (T  =  600  K)  where 

collisional  interaction  becomes  stronger,  the  E  model  cannot  be  fitted  to  the 

experimental  curve  of  Tr  =  600  K  for  any  reasonable  parameter  of  <AE>.  For 

each  model,  Fig.  4  shows  the  normalized  transition  probability  distribution 

P .  _  for  up-  and  down-transitions  from  the  level  energy  En  to  the  level  i  at 
1  L0  u 
a  representative  condition,  Tr,Tc  =  800,  400  for  the  parameters  obtained  by 

the  least  squares  fitting. 

It  is  of  interest  to  examine  the  relative  contribution  of  the  various 

elements  of  Nc  to  the  vibrational  activation  to  the  level  Eg  in  Nj .  The 

activation  distribution,  pF  ..jO vs_  E,  is  shown  for  the  temperature  combi na- 

L0  c 

tion  (800,400)  in  Fig.  5  for  the  various  model  parameters  fitted  to  experiment 

by  least  squares.  The  activation  functions  are  mainly  distributed  around  the 

initial  energy  2000  cm  .  The  maxima  would  be  shifted  to  higher  energies  if 

.  1 4 

<AE'>f  values  were  smaller  as  was  illustrated  previously  for  CB. 

L0 

The  various  transient  population  vectors  N-j  that  arise  after  single  colli¬ 
sion  are  shown  in  Fig.  6  for  the  E,  G  and  BE  models.  For  E  and  G,  the  popu¬ 
lation  elements  are  distributed  close  to  those  for  T^  =  800  K  in  the  energy 
4  -1 

range  below  E  =  10  cm  ,  while  those  for  the  BE  model  are  much  below  the  strong 
collider  curve.  The  relative  population  in  the  energy  range  below  Eg  provides 
the  main  criterion  for  the  average  amount  of  energy  transferred  on  collision, 
AEav  and  for  the  accommodation  coefficient  a  .  Thus,  the  BE  model  gives  rela¬ 
tively  small  values  of  AEflv  and  a  for  the  (Table  4).  The  same  situation  was 

14 

found  in  the  CB  system. 

By  contrast,  P values  are  determined  by  the  reacting  state  population 
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above  Eg,  in  the  vector,  and  which  derives  from  the  product  of  the 
up-transition  probability  to  states  above  Eg  with  the  state  populations  of 
the  Nc  vector  below  Eg.  And  although  a(BE)  is  smaller  than  a(G)  or  a(E),  it 
is  evident  from  Fig.  6  that  the  BE  model  is  relatively  very  efficient  in 
transferring  molecules  to  levels  above  Eg.  Indeed,  it  gives  rise  to  too-large 
values  of  Pc  at  higher  temperatures.  As  discussed  in  our  earlier  CB  study, 
the  models  for  P  that  fit  the  experimental  Pc  values  are  not  necessarily 
adequate  to  predict  a.  It  is  clearly  evident  that  measurements  of  a  and  of  Pc 
are  both  desirable  in  order  to  deduce  the  most  apt  form  of  P  . 

Plots  of  Pc(Tr,273)  obtained  by  simultaneous  measurements  of  gas  mixtures 

containing  both  1-MC8  and  CB  are  shown  in  Fig.  7.  Filled  squares  represent 

1 4 

the  P£  values  from  the  earlier  CB  experiment.  Good  reproducibility  for  CB  is 
evident  and  supports  the  presumption  of  a  negligible  role  of  catalytic  reaction 
on  the  fused  quartz  surface. 

Strong  collider  behavior  is  found  for  both  CB  and  1-MCB  molecules  curves 
at  Tr  <  450  K.  This  behavior  for  CB  has  been  shown  in  earlier  CB  studies. 

This  trend  is  related  to  the  fact  that  as  Tr  is  lowered,  the  time  that  molecules 
spend  trapped  in  the  gas-surface  potential  well  increases. 

Some  comparison  of  1-MCB  with  CB  for  the  single  collision  conditions  are 
given  in  Tab^e  6.  The  relatively  larger  experimental  values  of  Pc  for  1-MCB 
seems  to  result  from  the  difference  in  molecular  vibrational  eigenstate  densi¬ 
ties,  as  illustrated  by  the  mean  vibrational  energies  at  400  K:  1-MCB,  E4Qg  = 
1080  cm  CB,  E^gg  =  760  cm  ^ . 

The  least  squares  fit  values  of  <AE‘>F  for  the  G  model,  for  1-MCB  and  CB, 

fc0 

are  also  listed  in  Table  6. 

Q 

In  the  previous  VEM  studies  (m  <  1 )  of  Pc(Tr,273),  comparison  of  methyl- 
cyclopropane  with  cyclopropane  and  of ^ ^  cyclopropane  with  cyclopropane-dg  was 
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made  with  use  of  the  best  fitted  E  model  in  the  reactor  temperature  range  above 

T  =  800  K.  It  was  found  that  <AE‘>I-  decreased  with  an  increase  in  the  vibra- 
r  h0 

tional  heat  capacity  of  each  pair  of  molecules.  The  energy  transfer  (m  =  1 ) 
parameters  for  1-MCB  and  CB  fitted  to  the  experimental  pc(Tr>  2 73) values  are  listed 
with  use  of  the  G  model  in  Table  6.  The  more  complex  1-MCB  molecule  appears  to 
exhibit  somewhat  stronger  collider  behavior,  unlike  the  methyl-  and  the  cyclo¬ 
propane  systems  where  the  smaller  molecule  gave  larger  <AE'>  values.  The 

same  trend  is  observed  in  the  <AE'>C  values  calculated  with  use  of  the  best 

t0 

fit  E  model  to  the  Tp  =  800  curve, and  in  the  collisional  efficiency  .  We 
must  leave  this  effect  of  structure  for  future  resolution. 

The  G  model  leads  to  a  ^  0. 9-1.0  for  all  T^  used  (Table  4). 


Acknowledgment:  We  thank  Professor  H.  M.  Frey  for  the  gift  of  a  sample  of 
1-MCB.  This  work  was  supported  by  the  Office  of  Naval  Research. 
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Table  1.  Experimental  reaction  probability  per  collision  and 

collisional  efficiencies  in  1-methyl  cyclobutene  system. 


T\Tc 

273 

400 

500 

550 

800 

Pc 

2. 8x1  O'4 

3.  IxlO-4 

4.5xl0~4 

5.5xl04 

•l* 

0.035 

0.039 

0.056 

0.068 

700 

p 

c 

7.0xl0'5 

8.0xl0"5 

1.2xl0"4 

el 

0.081 

0.092 

0.138 

600 

P 

c 

6.5xl0'6 

8.6xl0'6 

1.25xl0"5 

el 

0.201 

0.265 

0.386 

a)  Collisional  efficiency  61  defined  as  6^  =  Pc(Tr,Tc)/P^c(Tr) , 

where  1 

s  c 

5C  (Tr)  is  the  reaction  probability  for  strong  collider 

interaction  of  the  molecule  with  the  surface. 


"vvc.  **>  • 
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Table  2.  Energy  transfer  parameters  fbrE,  G,  BE,  and  BG  models  obtained 
by  least  square  fitting  in  1-methyl  cyclobutene  system. 


Tr(K) 

Model 

<AE>  or  AE 

mp 

<AE 1 >r 

L0 

AEav 

800 

3280 

2850c 

890 

700 

E 

5820 

3650c 

800 

600 

___  b 

— 

— 

800 

3890 

4300C 

630 

700 

G 

4800 

5180C 

580 

600 

7220 

6430C 

580 

800 

1790 

3440 

1130 

700 

BE 

1730 

4780 

910 

600 

1770 

6850 

690 

800 

3570 

5210 

500 

700 

BG 

3570 

6290 

500 

600 

4160 

7780 

500 

The  energy  unit  is 

cm”1 

The  E  model  could 

not  be  fitted  to  the  600  K  curve 

in  Fig. 

Differ 

from  column 

3  due  to  truncation  effects 
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Table  3.  Average  energy  of  the  thermal  Boltzmann  population 
in  1 -methyl  cyclobutene  system. 


T(K) 

273 

400 

500 

600 

700 

800 

E(cm_1 ) 

410 

1080 

1870 

2860 

4030 

5320 

1 

14 

Table  4.  Average  energy  transfer  AEava  and  vibrational  energy  accommodation 
coefficient  ab  for  1-methyl  cyclobutene  by  single  collision  with 
a  seasoned  fused  quartz  surface. 


Tr(K) 

Model 

273 

400 

500 

550 

800 

4910,  0.94 

3630,  0.86 

2600,  0.75 

2090,  0.70 

700 

E 

3460,  0.96 

2470,  0.84 

1540,  0.72 

-- 

600 

-- 

— 

-- 

800 

4670,  0.95 

3950,  0.93 

3120,  0.91 

2640,  0.89 

700 

G 

3580,  0.99 

2870,  0.97 

2030,  0.94 

-- 

600 

2440,  1.0 

1720,  0.96 

910,  0.92 

— 

800 

990,  0.20 

960,  0.23 

860,  0.25 

780,  0.26 

700 

BE 

960,  0.27 

860,  0.29 

690,  0.32 

600 

920,  0.38 

730,  0.41 

440,  0.45 

a)  *Eav-Ef- 

Ec  (cm-1) 

cr 

P 

ii 

iEa/£r 

O 

1  LU 

1 

. tflCTTT 
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Table  5.  Some  comparisons  of  1 -methyl  cyclobutene  with  cyclobutene 
in  single  collision  systems. 


14 

CB 

l-MCB3 

Pc(800,  400) 

4xl0"5 

3xl0~4 

Pc(600,  400) 

lxlO"6 

9xl0"6 

Efl(kcal  mole"1) 

32.4 

34.2 

k(E)  (sec"1)5 

5 

10 

*  103 

fd 

1 

0. 64(800, 500)6 

E400^cm  ) 

760 

0.60(600,  273)C 

1080 

d 

<AE‘>f 

b0 

4090 

4300 

a)  Present  work 

b)  The  isomerization  rate  constant  by  RRKM  theory  at  several  100  cm 
above  Eq. 

c)  Average  thermal  energy  at  400  K. 

d)  Average  down  transition  energy  from  the  level  of  energy  EQ  , 
given  by  the  G  model  for  Tr  =  800  K. 


e)  Calculated  by  G  model. 
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a 

Table  6.  Energy  transfer  parameters  fitted  to  experimental  P  (T  ,273) 

C  t 

with  the  G  model  in  the  1-MCB  and  CB  systems. 


Tr(K) 

600 

700 

800 

<AE'>f 

t0 

CB 

5595 

4580 

4070 

(2490) b 

1-MCB 

6325 

5350 

4630 

( 2870) b 

AE 

a  v 

CB 

1620 

2420 

3200 

-O 

o 

CO 

CM 

1-MCB 

2440 

3590(3450)b 

4740 

( 4640 ) b 

CB 

0.96 

0.96 

0.94 

( 0 . 7  3) b 

a 

1-MCB 

1.0 

0.97(0.96)b 

0.97 

(0. 95)b 

CB 

0.13 

0.044 

0.021 

61 

1-MCB 

0.21 

0.084 

0.036 

a)  In  units 

of  cm”^ 

b)  Parenthetic  calculation  made  with  the  least  squares  fitted  E  model. 
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Fig.  1, 

Fig.  2. 
Fig.  3 

Fig.  4 

Fig.  5 


Figure  Captions 


Observed  isomerization  rate  constants  of  1 -methyl cyclobutene  to 

isoprene.  Open  and  solid  points  represent  the  total  rate  constants, 

k.,  and  the  wall  colli  si onal  rate  constants,  k  .  The  crossed  points 

were  not  used  for  the  P  calculations. 

c 

Plot  of  Pc  verus  the  initial  temperature  Tc- 

Calculated  P£ ( Tr ,Tc )  plots  for  each  finger  surface  temperature  by 

least  square  fitting  with  use  of  the  E  ( - ),  G  ( —  - ), 

BE  ( - ),  and  BG  ( - )  models;  the  heavy  solid  line 

summarizes  the  experimental  curve. 


.  Plot  of  normalized  transition  probability  distributions  for  down  and 
up  jumps  starting  from  the  vibrational  energy  level  EQ  =  11950  cm-1, 
calculated  for  the  E,  G,  BE,  and  BG  models  for  the  temperature  combi¬ 
nation  Tr,  Tc  =  800,  400. 


Curves  showing  relative  contribution  of  the  various  elements  of  Nc 

to  vibrational  activation  to  the  level  En  in  N, ,  given  by  Pr  .n  .. 

0  ~1  3  J  Eq-m  ci 

These  illustrative  plots  were  calculated  for  the  reaction  temperature 
combination  (800,400)  on  the  basis  of  the  least  squares  fitted 
parameters  and  are  plotted  with  arbitrary  units  of  the  ordinate. 
Obviously,  slightly  different  curves  would  result  from  a  different 
temperature  combination  or  different  fitting  procedure. 
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Fig.  6 


Fig.  7 


Relative  population  distribution  vector  versus  energy  E.  Nceci  and 
Nre<^  represent  the  Boltzmann  distribution  curves  characteristic  of 
the  initial  temperature,  Tc  =  400  K,  and  the  reactor  surface 
temperature,  Tr  =  800  K,  respectively.  N-j  superscripted  E,  6,  and 
BE  represent  the  transient  distribution  curve  after  one  collision 
with  the  surface,  the  two  left  ordinates  apply.  The  k(E)  curve  is 
the  RRKM  isomerization  rate  constant  to  which  the  right  side  ordinate 
applies. 

Plots  of  experimental  values  of  Pc( , 273)  versus  for  1-MCB  ( -  ) 

and  CB  ( - )  substrates  measured  by  internal  comparison  in  the  same 

conditional  reaction  system.  Also  shown  are  the  calculated  strong- 

collider  P^c (Tr )  curve.  Solid  squares  represent  earlier  experimental 
14 

values. 
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